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Abstract Wc investigate the effect of a variable, i. e. time-dependent, baekground on the standing acoustic 
(i.e. longitudinal) modes generated in a hot coronal loop. A theoretical model of ID geometry describing 
the coronal loop is applied. The background temperature is allowed to change as a function of time 
f\j I and undergoes an exponential decay with characteristic cooling times typical for coronal loops. The 

magnetic field is assumed to be uniform. Thermal conduction is the dominant mechanism of cooling the 
^5 I hot background plasma in the presence of an unspecified thermodynamic source that maintains the initial 

equilibrium. The influence of the rapidly cooling background plasma on the behaviour of standing acoustic 
(longitudinal) waves is investigated analytically. The temporally evolving dispersion relation and wave 
amplitude are derived by using the WKB theory. An analytic solution for the time-dependent amplitude 
that describes the influence of thermal conduction on the standing longitudinal (acoustic) wave is obtained 
f^ . by exploiting the properties of Sturm-Liouville problems. Next, numerical evaluations further illustrate 

CO I the behaviour of the standing acoustic waves in a system with variable, time dependent background. The 

results are applied to a number of detected loop oscillations. We find a remarkable agreement between 
^> I the theoretical predictions and the observations. The cooling of the background plasma due to thermal 

TY) • conduction is found to cause a strong damping for the slow standing magneto-acoustic waves in hot 

coronal loops in general. Further to this, the increase in the value of thermal conductivity leads to a 
strong decay in the amplitude of the longitudinal standing slow MHD waves. 



Q ' Keywords: Magnetohydrodynamics(MHD) • Plasmas • Sun: corona • Waves 



1. Introduction 

Recent consecutive solar observations by high-resolution imaging space telescopes and spectrometers have 
shown that the solar atmosphere is dynamic in nature and is composed of numerous magnetic structures of 
which coronal loops arc of the centre of focus in this paper. It has been confirmed that these complex struc- 
r — ' turcs of the solar corona can support a wide range of magnetohydrodynamic (MHD) waves which are nat- 

\^ . ural carriers of energy and may have the key to solve the problem of solar coronal heating (see, the recent 

-Y-\ ' review by e.g. Taroyan and Erdelyi) 120091 [McLaughlin et al. I201ip . In particular, one MHD wave mode 



f^ , present in coronal loops has become the centre of attention: the slow magneto-acoustic mode. The slow 

CN ' (propagating and standing) MHD waves are extensively supported by the solar coronal structures and are 

observed to be rapidly damped ( |Wang et a/., 2003a[ [Wang et al, 2003b[ |De Moortel, 2009[ [Wang, 2011] ) . 



Propagating intensity disturbances were first detected by the Ultraviolet Coronagraph Spectrometer 



onboard the Solar Heliospheric Observatory (SOHO/UVCS) along coronal plumes (Ofman et al., 1997 



Ofman et al., 1999 Ofman et al., 2000a[ ), and were identified as slow magneto-acoustic waves ( Ofman et al., 19^9| .[ 



C^ . Subsequently, similar intensity disturbances were observed in coronal loops by the Transition Region 
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and Coronal Explorer (TRACE) ( [Nightingale et al, 19991 |Schrijver et ai, 1999| |De Moortel et ai, 2000] 
|McEwan and De Moortel, 2006[ ) and the EUV Imaging Telescope onboard the Solar Heliospheric Obser- 
vatory (SOHO/EIT) dBerghmans and Clette, 1999[ ). iNakariakov etaD ([20001 found that the slow MHD 
wave evolution is affected by dissipation and gravitational stratification. [Erdelyi and Taroyan| (j2008p and 
[Wang et al.\ (j2009p have detected longitudinally propagating slow MHD waves with around five-minute 
period in the transition region and five coronal lines at the footpoint of a coronal loop by Hinode/ElS. 
The source of these oscillations are suggested to be the leakage of the p modes from photosphere region 



through the chromosphere and transition region into the corona ( De Pontieu et ai, 2005 1 . 

Oscillations interpreted as longitudinal standing (slow) magneto-acoustic waves have been observed 
in hot (T > 6 MK) active region loops by the Solar Ultraviolet Measurement of Emitted Radiation 



(SUMER) spectrometer on board SOHO ( |Wang et ai, 2002| [Wang et ai, 2003b[ [Taroyan et al, 2007D . 



These oscillations have periods in the range of 8.6 to 32.3 minutes with decay times of 3.1 — 42.3 
minutes and amplitudes between 12 and 353 km s~^ ( Wang et ai, 200"5| . Moreover, IMariskal (j2005p 



have reported Doppler shift oscillations during solar flares with Yohkoh in a high temperature region 
reaching 12 — 14 MK. These oscillations are interpreted in terms of the standing-slow mode MHD waves 
( Mariska, 2006D . Evidence for the standing slow mode can be underpinned from the phase relationship 



between velocity and intensity where a quarter-period phase difference is a characteristic of the standing 
waves while the propagating waves exhibit an in-phase relationship. Therefore, an approximate quarter- 
period delay of the intensity variations behind the Doppler shift strongly support that the oscillations 
observed by SUMER are slow standing modes. The observed oscillations in coronal loops indicate that 
the standing slow modes are likely triggered by micro-flares which are produced by impulsive heating 
( Mendoza-Briceho et ai, 2002p. In a recent work,[Taroyan and Bradshaw (|2008p found that longitudinal 



standing waves can be excited in cooler EUV loops under the effect of all important mechanisms such as 
gravitational and thermal stratification, losses, etc. According to this study, hot loops are not only the 
origin of the standing waves but these waves can also be formed in cooler loops. 

Nowadays, the damping of slow magneto- acoustic waves have become a subject of remarkable obser- 
vational and theoretical attention due to a possibility of damping timescalc of the waves on disclosing 
the physical processes that dominate the energy of the coronal loop in which they are detected. The 
majority of theoretical and numerical studies on damping of propagating and standing slow MHD waves 
show that the understanding of dominant mechanisms of rapid damping can be captured from ID 
linear ( Sigalotti et ai, 2007P or nonlinear model ( [Wang, 201ip . For instance, |Ofman and Wang] (|2002p 



and IMendoza-Briceho et al.\ (|2004p found that the standing MHD waves are strongly damped because of 
thermal conduction in nonlinear model whereas using a linear MHD model [Pandey and Dwivedi| (|2006p 
indicated that the individual influence of thermal conduction or viscosity is not enough to account for 
the observed damping. In a static ID isothermal medium, IDe Moortel and HoodI (|2003p investigated the 
behaviour of both propagating and standing slow MHD waves in the presence of thermal conduction 
and compressive viscosity. They found that thermal conduction is the dominant damping mechanism of 
thermal and magneto-acoustic waves in coronal loops. It is also shown that thermal mode is damped in 
the form of standing waves. Moreover. IMorton et aLI(|2010p have shown that the radiative cooling causes 
a damping of the slow mode by up to 60% within characteristic lifetimes. 

IMendoza-Briceho et aLI(|2004p studied the influence of gravitational stratification on damping of stand- 
ing MHD waves in hot coronal loops and found that enhanced nonlinear viscous dissipation due to gravity 
may reduce the damping times by about 10 — 20% compared to the unstratificd loops. In contrast to 
most recent work, [Sigalotti e^"ar](|2007P found that thermal conduction can only be accounted as damping 
mechanism when the compressive viscosity is added to the model. [Bradshaw and Erdelyi[ (j2008p reported 
that the radiation due to a non-equilibrium ionisation balance could cause by up to 10% reduction of 
wave-damping timescale in comparison to the equilibrium case. IVerwichte et al.\ (j2008p showed that shock 
dissipation at large amplitudes gives rise to enhancement of the damping rate which is up to 50% larger 
than given by thermal conduction alone. In non-isothermal, hot, gravitationally stratified coronal loops, 
[Erdelyi et al.\ (j2008p investigated the damping of standing slow (longitudinal) waves and found that the 
decay time of waves decreases with the increase of the initial temperature. Further to this, they derived 
a relation between the damping time and the parameter determining the apex temperature in the form 
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of second-order scaling polynomial. In all these earlier studies the background (equilibrium) plasma was 
considered time-independent. This is certainly a working model as long as the background plasma changes 
occur on timescales much longer than typical perturbation timescales. In a recent theoretical work by 
[Erdelyi et aL|(j201ip . the cooling of the background plasma due to thermal conduction has been found to 
decrease the rate of damping of the propagating magneto-acoustic waves in a weakly stratified atmosphere 
where temporal changes in the slowly evolving equilibrium were considered. 

In this study, we examine how the cooling background state affects the longitudinal standing waves in 
a uniform magnetised plasma. The additional complexity of a time-dependent background gives us a step 
along to model hot coronal loops. The cooling of the background plasma is assumed to be dominated 
by thermal conduction in the presence of an unspecified thermodynamic source where the plasma is 
assumed to be cooling in time with an exponential profile. The coefficient of thermal conductivity is 
presumed to be varying as a function of time to investigate how the variation of thermal conduction 
influences the behaviour of standing slow-mode waves. The geometry of coronal loop is described by a 
semi-circular shape. An analytical solution for the time-dependant amplitude of the standing slow MHD 
waves is derived using the WKB theory with the aid of properties of Sturm-Liouville problems. The 
results exhibit that thermal conduction causes a significant variation on the amplitude of the standing 
magneto-acoustic modes with time. The amplitude of waves are found to suffer damping due to the 
cooling of the background plasma by thermal conduction. The influence of the ratio change of the period 
of standing oscillation to the cooling time scale on the rate of damping of hot-loop oscillations seems to 
be highly efficient, where increasing the value of this ratio deceases the strength of damping rapidly. 

2. The Model and Governing Equations 

Consider a uniformly magnetised plasma in which the temperature is changing as a function of time due 
to thermal conduction and the density is a constant. The magnetic field is assumed to be uniform and 
aligned the z-axes, i.e. Bq = Bqz. Therefore, the background state can be described as follows: 
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Figure 1.: Coronal loop 
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Here Tq, po, po and Bq are the background quantities identifying the temperature, pressure, density and 
magnetic field, respectively; and P is the period of the loop oscillation and Tc is the cooling time scale. 



The governing MHD equations for the background plasma take the following form 

|+V.(pv).0, 

dv 1 

p—- + p(v.V)v = -Wp + —(V X B) X B, 
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dt 
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where v is the flow velocity, B is the magnetic field, g is the gravity, /io is the magnetic permeability 
of free space, 7 is the ratio of specific heats, R is the gas constant, jl is the mean molecular weight, T 
is the temperature, V(«;||VT) is the thermal conduction term where K|| = kqT^/^, C is an unspecified 
thermodynamic source term, p and p are the plasma density and its pressure, respectively. 

Assuming that there is no background flow and the background density is constant, so the previous 
equations determining the background plasma state reduce to 



uo = 0, Po 
Vpo = 0, 

R 

Po = —poTo, 
P 

R dTo ^ 
jj, at 



const. 



i.e.po ^Tq, 



(11) 
(12) 

(13) 
(14) 



where the index denotes background quantity. The observed cooling of coronal loops has been shown 
to be well-approximated by exponential profile for radiative-cooling loops (jAschwanden and Terradasl 
[20051 [Morton and Erdelyi| [^OOUl [20TU| . More recently, [Erdelyi et d] (P0TT|) found that the temperature 
decreases exponentially in time due to thermal conduction in hot coronal loops. As a result, we will 
assume here that the cooling temperature profile of coronal loops has the form 



To = ro,exp( ), (15) 

Tc 

where Tgi is the initial temperature at t = 0. 

Let us perturb the background state, where all variables can be written in the form 

/(z,i) = /o(t) + /i(z,t), 

where the subscript refers to the equilibrium quantities and the subscript 1 indicates the perturbed 
quantities. In this study we consider longitudinal motions only, i.e. v — vz, so the linearised perturbed 
MHD equations in a ID system are 



(16) 
(17) 
(18) 
(19) 
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Here vi = viz is the longitudinal velocity perturbation. It is clear that there are no terms exhibiting the 
magnetic field but the standing waves are still guided by the magnetic field. Non-dimensionalisation will 
be used to simplify the ID governing Equations (J16p — (|19p . The dimensionlcss quantities are introduced 
as 
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(20) 



where the subscript i represents the value of the quantity at t = 0, L is the loop length, and Cgi is the 
initial sound speed. 
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Now, we aim to find the governing equation for tlie perturbed velocity, wliicli leads to reveal the 
behaviour of the standing magneto-acoustic mode subject to initial conditions generated in a hot coronal 
loop. Using the continuity and ideal gas law equations, Equation (|18|) in terms of dimensionless variables, 
removing the tilde for the sake of simplicity, takes the following form 






Tt 



d'^Vl — (7 



dz^ 



= T, 
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(7 - l)/^^o ^oY^ 
RLy/jpoipo 



(21) 



where a is defined as thermal ratio and found to be a small quantity under typical coronal conditions 
(see, e.g. IDe Moortel and Hoodi I2003P demonstrated here as 



f To = 6 X 10^ - 5 X 10^ K, 
po = 1-67 X 10-12 kg m-^ 
K|| = 10-" To^^ W m-i deg-\ 

<; /i = 0.6, 

i? = 8.3x lO^m^ s-2 deg-\ 
7 = 5/3, 

\, L^IO^ m, 

giving characteristic values of cr € [0.0068, 0.48] for T € [6 x 10^ 5 x 10^] in K. 
Equation dH]) with the aid of Equations (HI]), (US]) and ([H]) will be 
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Differentiating with respect to time and, using Equations p6[) and (|17p . we obtain 



(22) 
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which is the governing equation and can be recast to a convenient form for further analysis: 

^3/2 QJo o^vi 
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(25) 



There are three different cases with interest for the present context one can recover from Equation (^51 



Case I. In the absence of the thermal conduction, a, and the unspecified thermodynamic source in the 
energy equation, C oc dTo/di, the governing equation reduces to 



d'^Vi 2^^'^l 



dt^ 



dz'^ 



0, 



'To = const. 



which has the solution 



vi{z,t) = acos(7rz) cos(7rCsi), 



(26) 



(27) 



under the boundary conditions representing a line-tied fiux tube perturbed as a semi-circular initial 
motion proportional to its fundamental mode 



1 ^ 

ui(t,±-) = 0, wi(0, z) = acos(7rz), — !-(0,z) = 0. 



(28) 



where a is the initial amplitude of the standing wave at t = 0. A more general, e.g. broad-band, 
perturbation would, of course, give the solution in the mathematical form of a Fourier-series. 
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Case II. In the case of no thermal conduction, i.e. a = 0, the effect of unspecified thermodynamic source 
in the system will be found by solving the following equation 

|2„ 



d Vl 2<^ ^^i 



Qt2 ^'^=0' Cs{t) = ^%{£) ^ const., (29) 

which is formally exactly Equation (j26|) but with variable background temperature, Tg ~ T(){t). In 
spite of the absence of a, the coefficient of the bracket in the first term in the right-hand-side of 
Equation (|25p is originally derived from thermal conduction term as seen in Equation p3p . Therefore, 
this term is not added to Equation (P^ . 

Case III. The effect of thermal conduction on the behaviour of the standing acoustic waves in the 
presence of an unspecified thermodynamic source will be investigated by solving Equation (|25p . 

3. Analytical Solutions 

Our goal now is to find an analytic solution to the governing Equation (j29p first in case II and next to 
Equation ([25]) in case III. Let us point out that Equations (|25p and ([29)) have derivatives multiplied by 
small factors a and e so this enables the use of the WKB theory to obtain an approximate solution, where 
the WKB estimates give good approximations when the smaller value of factor used. 

Case II. Let us first start to solve Equation (|29p to establish the behaviour of standing magneto- 
acoustic waves under the influence of an unspecified thermodynamic source in a cooling background state. 
Assuming that ti = et which is defined as a slow timescale, meaning the cooling timescale is (much) longer 
than the period of the oscillations. Equation ([^ will reduce to 

In line with applying the WKB estimates let the perturbed velocity variable have the form 

?;i(2,ii)=Q(z,ii)expQe(ti)j . (31) 

The amplitude Q can be expanded in the power series as follows 

Q(z,ii) = Oo + eQi + ---. (32) 

Substituting Equations ((3T|) and (|32p into Equation (pO)) . and taking terms of order e^^ we obtain at 
leading order 

^ + ^^0=0, (33) 

where w = dO/dii. The boundary conditions (j28p applicable to the function Qo is 

Qo = at ^ = ±^- (34) 

Equations (|33p and p4p represent a boundary-value problem that determines the frequency of the 
standing longitudinal (acoustic) mode in a cooling plasma with a varying temperature as function of 
time. The general solution physically acceptable to this boundary value problem has the form 

oo 

Qo(2,ii) = ^^n(ti)cos((27i-|-l)7rz), a;„ = (2n + l)7rcs, 7i = 0, 1, 2, • • • . (35) 

n=0 
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Then, collecting terms of order e^^, we obtain the equation determining the amplitude 






n = ^ 



—- Qo + 2a;—— 
oil ati 



(36) 



It follows from Equation 



that Qi satisfies the boundary conditions 
Qi=0 



at z 



4 



(37) 



The boundary- value problem, Equations p6p and (|37p . determining Qi is a Sturm-Liouville problem and 
has a solution only when the right-hand-side of Equation p6p satisfies the compatibility condition, which 
is the condition that it is orthogonal to Qo (see, IRudermaiil I201ip . This condition can be obtained by 
multiplying Equation p6p by Qo, integrating with respect to z from —1/2 to 1/2 and using the boundary 
conditions (j37p . The compatibility condition is eventually written as 



1/2 



1/2 



dii oil 



dz = 0. 



(38) 



which gives the amplitude of the standing wave in the following form 



,ii 



c.(0) 



An (ii ) - A„ (0) exp( I ) - A„ (0) ^ ^^ ^^^ ^ 



(39) 



The value of constant A„(0) can be found from Equation (|55|) at ti = 0, which is in the form of Fourier 
cosine series, using the boundary condition (PSJ) . Then, the solution ([M)) . in scaled (i.e. physical) variables, 
takes the form 



An{t) = a. 



^ c,(0) 



(40) 



This result describes the variation of a time-dependent amplitude of longitudinal standing waves in a 
dynamically cooling magnetic flux tube. Apparently, Equation PO)) indicates that the cooling (or heating) 
leads to an amplification (decrease) of loop oscillation. 



Case III. Next, the behaviour of the standing wave in a system dominated by thermal conduction 
will be obtained by solving the governing Equation ([^5)) . Similarly, the WKB theory will be used to find 
the solution of Equation (1251) . Let us now introduce two slow variables ti ~ et and a — ea, so Equation 
(PSj) will be 



Q3 



Vl 



Id'^vi 



53, 



Vl 



dtl 2 dtl £2 Qi-^Q.^2 



~ 5 d^vi 
(Tc: 



d'^vi 5 c^ d^vi a c^ d'^v\ 



• dhdz^ 



' dtjdz^ 2 £2 dz^ £2 7 dz'^ 



0, 



(41) 



and the perturbed velocity by the new scaled variables and the WKB approximation is given by Equation 
([5T|) . Substituting Equations ([31]) and ((5^ into Equation (PTj) . and taking the highest order terms in e, 
which is again £~^, we obtain 

(42) 



d^Qo 



9z2 






with the boundary conditions 



Qo = 



at z = ± — . 
2 



(43) 
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The solution to this boundary-value problem is given by 

oo 

Qo{z,h) = Y^Bn{ti)cos{{2n + l)TTz), a;„ = (27i + l)7rc,, n = 0, 1, 2, • • • , (44) 

n=0 

where Bn{ti) stands for the amplitude of the longitudinal standing modes and it will be found by taking 
the second highest order terms in e for Equation (j41[) and then using the properties of Sturm-Liouvillc 
problems as follows. 

The next largest order terms in e, of order e~^, give the following equation 

d'^Qi , uj'^ ^ i \,7 2 , o dw j^Qo , 2 9^Qo , ,5 2 ~ 2 5n^^Qo ~_cl d'^Qo 



dz^ ci u)6 



i-u + 3uj-—)Qo + 3uj —-+ c^^r—-^ + i-c: - ^uj cj— -^ - a- 



-2 dt/^^ dh '"dhdz^ "■2'" "' dz^ 7 dz^ 

(45) 
It yields from Equation ((28)) that Qi satisfies the boundary conditions 

Qi=0 at 2 = ±^- (46) 

Analogous to Equations ([55)) and (|57)). this boundary- value problem, Equations (|15)) and (HHj), has a so- 
lution only when the right-hand side of Equation (|45p satisfies the compatibility condition. Consequently, 

-1/2 ^C2 



^'''2,0 '^^ \/^'2 I Q 2^ ^Qo , 2^ d^Qo I /5 2 ~ 2 5\^ d'^Qo 



r.7 



c: „ d^Qa 



^-Qo- 



7 ^" 9z4 
Substituting (|44| into ()47[) , we obtain the amplitude of standing wave 

1. , ^^7-l^,o„ , -,^2_2/„5 



dz = 0. (47) 



B„ih) = B„(0)exp ( -ti + -(^)(2n -H l)'^'(c^ - 1)1 , (48) 

which can be re-written in the scaled variables as 

Bn{t) = i?„(0) exp f^t + ^(l^)(2,i + 1)2^2(^5 _ ^^ (49) 

Now, Equation (j^l)) . which represents a Fourier cosine series, with the boundary conditions (|^5)) will be 
applied to obtain the value of constants -B,i(0). Hence, the solution ([^^ is 

i?„(i) = aexp (U + f (l^)(2n + lfn^{4 - 1)) , (50) 

\4 5e 7 / 

which reveals the temporal evolution of longitudinal standing-mode amplitude due to thermal conduction 
in a temporally variable (cooling or heating) background plasma. Note that, in the limit cr — >■ 0, i.e. if 
there is no thermal conduction. Equation ()50p reduces to Equation ()40p which represents the amplitude 
variation of standing waves generated in a inodel of merely cooling plasma without non-ideal effects as, 
e.g., thermal conduction. 



4. Numerical Evaluations 

IMorton et al.\ (I2010p and jErdelyi et al.\ (|201ip have shown that the WKB approximation estimates accu- 
rately the solutions to the frequency and amplitude variations in time and space for waves supported by 
oscillating magnetic loops due to plasma cooling by radiation and/or thermal conduction, respectively. 
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The obtained approximate solutions by the WKB theory can be motivated using numerical evaluations 
to demonstrate a clear view of the behaviour of MHD waves which is analytically found. 

The amplitude of the longitudinal (acoustic) standing waves is plotted after calculating the variables 
using standard coronal values. Figure [2] shows the variations of the amplitude due to cooling (or heating) 
mechanism. This graph exhibits an amplification for the wave amplitude caused by the cooling of the 
background where increasing e (the ratio of the period of oscillation to the cooling time scale) increases 
the amplitude. It should be mentioned that this result is in agreement with that reached by |Erdelyi et a/.] 
(|201ip who found that the efficiency of damping is reduced by the cooling background plasma, where 
Figure [5] displays no damping for standing oscillation because the model is purely dominated by the 
cooling. 
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Figure 2. The amplitude of the standing wave with different values of e (0.01, 0.05, 0.1, 0.5) representing the ratio of period 
to the cooling time. 



Next, in Figure |3] we demonstrate how the amplitude of standing longitudinal (acoustic) waves for a 
range of values of e and as function of the value of thermal ratio a is changing in different temperature 
regions. It is found that the variation of e leads to a considerable change in the rate of damping of 
both cool and hot loops. Figure I5al illustrates the trend of amplitude of the EUV (cool) loops such that 
the oscillation amplitude decreases slowly in regions of temperature 600, ODD K. The decay of the wave 
amplitude increases slightly with time in hot corona, for instance loops of temperature 3 MK and 5 MK 
as depicted in Figure [3bl and [5cl respectively. The strength of damping of hot {e.g. SXT/XRT) loops is 
found to be much stronger and shows little change for the smallest values of e, i.e. 0.01 — 0.1, whereas 
large enough values of e, i.e. in the range 0.1 < e < 0.5, cause a rapid reduction in the rate of damping 
of the wave amplitude. Moreover, it is obvious from Equation (|49p that the first term in the exponential 
function is the same as in Equation (|39p were dominated by the cooling. This indicates that the emergence 
of cooling by thermal conduction in the system of hot coronal loops decreases the rate of damping as 
exhibited in Figure [31 

In Figure 3] we present the effect of varying the magnitude of thermal conduction coefficient, k, on 
the rate of damping of the standing acoustic wave. Typical values for the coefficient of the thermal 
conductivity k = [10~^°, 10~^^, 10~^^]r^/^ are taken to shed light on the influence of thermal conduction 
on hot corona, T = 3 MK, and the characteristic value of e = 0.1 is assumed (see, IPriestl I2000p . It is 
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Figure 3. The amplitude of the standing wave with different values of e (0.01, 0.05, 0.1, 0.5) representing the ratio of period 
to the cooling time and specific value of u. i.e. the value of thermal ratio, (a) a = 0.0068 (T = 600000 K), (b) cr = 0.17 
(T = 3 MK), (c) o- = 0.48 (T = 5 MK). 

found that the rate of damping is changing rather rapidly with ahering k by just an order of magnitude 
where increasing the value of k gives rise to a strong decline in the amplitude of the standing slow mode. 
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Figure 4. The amplitude of the standing wave with different values of the thermal-conduction coefficient. 



K0=(10- 



,10" 



,10 ) and specific value of the ratio of period to the cooling time, e = 0.1 where T = 3 MK. 



The significance of the obtained results is to be comparable to observations. We only present here a 



quantitative comparison with the properties of hot loop oscillations observed by SUMER (see, Wang et al 
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I2003ap . where the standing slow- mode waves are detected only in the region of temperature > 6 MK. 
The periods of oscillations are 7 — 31 minutes. The slow standing wave has seen to be strongly damped 
with characteristic decay times 5.7 — 36.8 minutes mainly likely due to thermal conduction. The typical 
length of coronal loops is around 230 Mm. 

In our work, we found that hot loop oscillations experience a strong damping due to thermal conduction 
which might be comparable with the observed damping provided the value of e is small enough as shown 
in Figure [5cl Further to this, Figure |3] exhibits that the large value of thermal conduction coefficient leads 
to a rapid damping which is likely applicable for the observed damping of standing acoustic modes as 
discussed by |Ofman and Wang| ([^00^ and IMendoza-Briceno et al.\^^UU^ . 

5. Discussion and Conclusion 

In this work, we have investigated the influence of a cooling background on the standing magneto- 
acoustic waves generated in a uniformly magnetised plasma. Thermal conduction is assumed to be the 
dominant mechanism of cooling background plasma. The background temperature is allowed to change 
as a function of time and to decay exponentially with characteristic cooling times typical for coronal 
loops. The magnetic field is assumed to be constant and in the z direction which may be a suitable 
model for loops with large aspect ratio. A theoretical model of ID geometry describing the coronal loop 
is applied. A time-dependant governing equation is derived by perturbing the background plasma on a 
time-scale greater than the period of the oscillation. Three different cases were considered: (I) absence 
of thermal conduction a and unspecified cooling or heating mechanism £. (II) presence of unknown 
thermodynamic source C only, i.e. cr = 0, (III) the influence of thermal conduction a combined with the 
unknown mechanism C. 

The WKB theory is used to find the analytical solution of the governing equation in case II and III 
where the governing equation in case I is solved by a direct method and gives the undamped standing 
wave, i.e. sound speed is constant. An approximate solution that describes a time-dependant amplitude 
of the standing acoustic mode is obtained with the aid of the properties of a Sturm-Liouville problem. 
The analytically derived solutions are exhibited numerically to give much illustration to the behaviour of 
MHD slow waves. 

In the second case, the individual influence of cooling background plasma on hot-loop oscillation is 
found to cause an amplification to the amplitude of the longitudinal standing wave. It is noted that the 
rate of amplification varies according to the change in the ratio of the oscillatory period to the cooling 
time scale e, where increasing e increases the amplitude of MHD wave. 

In the third case, which is of our interest, the temporally dependant amplitude is found to undergo a 
strong damping due to the cooling of the background plasma by thermal conduction in the hot corona. 
Further to this, we note that the presence of cooling in a model of hot coronal loop decreases the efficiency 
of damping. The variation of the ratio of the period of oscillation to the cooling time scale, e, plays an 
effective role on changing the rate of damping of oscillating hot coronal loops, causing a fast decline in 
the decay degree of oscillation amplitude once the value of this ratio is so large. 

The obtained results indicate that this investigation contributes with the previous studies ( [Morton and ErdelyiH 
[^OM \SUW\ IMorton et 'aL[ [20111 [Erdelyi et aL\ [^UTT]) on demonstrating that the temporal evolution of 
coronal plasma due to the dissipative process, i.e. the cooling of the background plasma due to radia- 
tion/thermal conduction, has a great influence on the coronal oscillations. In the modelling of solar coronal 
loop, the temporal and spatial dependant dynamic background plasma is necessary to be considered to 
understand the properties of observed MHD waves. 
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